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Abstract 

We consider the specific models of Zliu-Kroemer and BenDaniel-Duke in a secli^- 
mass background and point out interesting correspondences witli tlie stationary 1- 
soliton and 2-soliton solutions of the KdV equation in a supersymmetric framework. 
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In dealing with position dependent mass (PDM) models controlled by a sech^- 
mass profile, we demonstrated [1] recently that, in the framework of a first-order 
intertwining relationship, such a mass environment generates an infinite sequence 
of bound states for the conventional free-particle problem. Noting that the inter- 
twining relationships are naturally embedded in the formalism [2] of the so-called 
supersymmetric quantum mechanics (SUSYQM), we feel tempted to dig this issue 
a little deeper by choosing to examine the connections between the discrete eigen- 
values of such a PDM quantum Hamiltonian (transformed appropriately so that a 
SUSY structure is evident) and the stationary soliton solutions of the Korteweg de 
Vries (KdV) equation that match with the mass function upto a constant of pro- 
portionality. 

Let us begin with the standard time-independent representation of the PDM 
Schrodinger equation [3] 

where M[x) is the dimensionless equivalence of the mass function m(x) defined by 
m{x) = moM(x), and we have chosen units such that h = 2mQ = 1. The effective 
potential Vg// contains, apart from the given V{x), the real ambiguity parameters 
a and /? whose occurrences are typical in PDM settings : 

1 M" M'^ 

Veff = V{x) + + 1)— - {a{a + p + 1) + p + 1}^ (2) 

Suitable physical choices of a and /? have been reported in the literature [4-25] but of 
particular interest to us are the schemes of Zhu-Kroemer (ZK) [4] [a = —1/2,(3 = 0] 
and BenDaniel-Duke (BDD) [5] [a = 0, P = —1] which were shown [1] to be dual of 
each other for the free-particle case V{x) = Vq that is independent of any choice of 
M{x). 

Substituting (2) in (1) and assuming for Vq the form 

T/q = e - A(A + A,geR (3) 

we can recast (1) to the standard constant-mass Schrodinger equation namely 

(-^ + «)V^ = (4) 
with the energy level term missing. In (4), u is given by 
-3 . . . . . ,iM'2 l^M" 



u 



- - {a{a + P + + + - A(A + l)Mq'. (5) 
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However, Eq. (4) can also be regarded as the linearized partner of the Riccati 
equation 

■u = f ^ + f ' (6) 
upon putting v — ^. The latter is the Cole-Hopf transformation. 

A nonlinear connection such as the one given by (6), also known as the Miura 
map, has an interesting implication. It transfers a solution of the modified KdV 
equation 

Vt = 6v\' - v'" (7) 
into a solution of the KdV equation 

Ut = 6uu' - u'" (8) 

which is straightforward to check. 

The KdV equation has a very rich internal structure [26-27]. In particular, it 
admits of a Lax representation Lt = [B, L], where L = —d^ + u is a Schrodinger-like 
operator and B is given hj B = —49^ + 6udx + 3u'. One can solve for L in the 
from L{t) = S{t)L{0)S^^{t) with St = BS. The related eigenvalue problem then 
implies that the spectrum of L is conserved and yields for the KdV an infinite chain 
of conserved charges. 



Noting that the KdV is invariant under the set of transformations 

t^t', x^x'-6ct\ u^u' + c (9) 
where c is a constant, the energy levels /in can be introduced in (4): 

{-^+U)ij = ^Ini^. (10) 

The manner of interplay between the PDM form (5) of u for specific choices of the 
parameters a, j3 and the initial condition u{x^ 0) = Uq used as inputs to solve for the 
KdV (as is normally done in the inverse scattering problem) is our point of enquiry. 



It can be proved that the discrete eigenvalues //„ are time-independent. For this 
we have to express the KdV in the conserved form 

Ut (-3it^ u^^)^ = (11) 

and substitute u from (10) into it. We obtain 

{y^nW + {i^(t>x - i^x4>)x = (12) 

where (p — ipt + i^xxx ~ 3(w -|- n)'il)x- On integrating (12) we find {iin)t — where 
we have employed normalized ^0 and considered vanishing asymptotic conditions for 
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■0 and its derivatives. The eigenvalues are determined using for the potential 
the initial value Uq that corresponds to a stationary soliton solution of the the KdV 
equation. 

In the context of (10), the Riccati equation (6) is transformed to 

u^v^ + v' + n (13) 
where v as a solution of the generalized MKdV equation 

= 6{v^ + ij)v' - v'" (14) 
ensures that u evolves according to the KdV equation. 

For the 1-soliton and the 2-soliton solutions of the KdV, the corresponding start- 
ing solutions Uq along with if], v and the eigenvalues fin (n=l, 2) are given by : 

1- soliton : x^q^^ = —2q^sech^qx, 

■01 = -^sechqx, v^^^ — —qteaihqx, jii — —q^ (15) 
v2 

2- soliton : Uq^ — — 6g^sech^qx , 

'02"'* — -^sech^qx, ^2"'* = — 2gtanhga;, 112'^ — —4q^ 



and 



(6) (b) 1-2 tanh^ qx (6) 2 

■02 = — secnqx ta.nhqx, V2 = q , l^, — ~<1 

2 tanh qx 



where note that for one discrete value of the Schrodinger equation (10), there exists 
a 1-sohton solution and vice- versa. Similarly for the 2-soliton case. Here the '0's are 
normalized. 

The results in (15), which can also be extended to the N-soliton case, have been 
obtained by solving the eigenvalue problem for the Schrodinger equation (10). The 
solutions Uq^ and Uq^ act in (10) as the refiectionless potentials. The inverse scat- 
tering method, which exploits this refiectionless feature, determines the evolution 
of the scattering parameters. Subsequently the Gelfand-Levitan integral equation is 
solved to obtain the solution u{x, t) of the KdV equation. 

Turning now to the PDM induced u given by (5), we immediately recognize from 
(10) that for the choice of the mass function M{x) — sech^qx, the ZK scheme yields 
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the 1-sohton results U(^\'ijji(iJ, — —q^) corresponding to A = 1, —2 and the 2-soliton 
result u'q\ ijj2'\lJ' — — 4g^) and ■02^^ (a* — ~(f') corresponding to A = 2, —3. 

On the other hand, the BDD scheme is consistent with the form 

w = ?2(l-2sech2?x), A =1,-2 (16) 

for -01 [ji — 0) and 

M = g2(i_6sech2ga;), A = 2, -3 (17) 

for both the sets '4>2'\l^ — "3^^) and '4>2\l^ — 0)- 

To interpret the above results, a few remarks on SUSY are in order [28]. We 
first of all verify that not only (13) but also u = v"^ — v' + ^ carries a solution of the 
generalized KdV (14) into a solution of the KdV. 

Denoting 

wc notice that the combination V^'^^ can be identified as the usual partner potentials 
of SUSYQM. 

To examine the role of V"*^'*'-' in the present context, let there be a Hamiltonian 
Hi with potential Vi that is asymptotically vanishing and having a set of n discrete 

eigenvalues i^i, IJ>2, , A*n- If we define = V\— fin+i then, in unbroken SUSY, 

we at once know that the spectra of V^~^^ and those of V'^^^ are one to one except 
that the latter has an additional /i = state. In other words, the eigenvalues of V^~^ 

are /xi — /^n+i; A*2 — l^n+i, , l^n ~ and 0. This means that the Hamiltonian 

H2 with potential V2 defined by V2 — V'^'^ + jin+i has (n+1) discrete eigenvalues 

A*l, A*2, , A*n, A'n+l- 

Let us apply the above ideas to the simple case of Vi = and generate the 
corresponding potential V2 with a single bound state with 11 = —q^ [29]. We have 
— v'^ — v' — q^ > in other words, has no bound state at all. Solving 
we get V = —qtanh qx (i.e. the 1-soliton result) which in turn gives V^~^ — q^{l- 
2sech^ga;) that supports a zero energy {/i — 0) bound state ■0o ~ sechqx : 

H_^ljo = iJo + {v'' + v')ilJo (19) 
Thus V2 — — 2g^sech^ga; has a single bound state. 

We immediately recognize V2 and V^^^ to be the PDM potential u for the ZK 
and BDD schemes respectively corresponding to the 1-soliton case. The same is 
true for the 2-soliton results with v matching with the 2-soliton solutions and 
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emerging similar to (17). 



One-dimensional supcrsymmetric approach to PDM quantum systems has been 
explored before in PDM scenarios. The partner potentails were found to obey [22] 
the same PDM dependence but in different potentials. The approach of this work 
is however different in spirit from such a viewpoint in that we have sought to es- 
tablish a link between a hierarchy of refiectionless potentials ( corresponding to the 
stationary soliton solutions of the KdV ) with an arbitrary bound state spectrum 
and those of SUSY in PDM models for suitable values of the ambiguity parameters. 
Our starting potential pertaining to the free-particle case V{x) — Vq can be made 
to coincide with V^'^^ by choosing, for example, e = 3q'^ in the 1-sohton case and 
e = 12q'^ in the 2-soliton case. 

Finally, we can extend our treatment to other special cases of the effective po- 
tential Veff namely those of the Bastard [30] and Li and Kuhn (redistributed) 
[31] Hamiltonians. For the 1-soliton result of (15), u for the Bastard scheme is 
u — —q^il + 3sech^ga;) (/i = — 2g^) while for the 2-soliton results given by (16), 
u turns out to be —q'^{l + 6sech^ga;) both for t/^""* '^2 \ with an associated ji- 
value oi n = —bq^ and /x = —2q^ respectively. However, in the Bastard model A is 
non-integral. A non-integral A also emerges in the Li-Kuhn scheme where we find 
u — — 2g^sech^5x(// = —q^) for the 1-soliton solution and u — — Gg^sech^^x for both 
the 2-soliton solutions of ■02"^ (a* — and '4'2\l^ — ~<f)- 

I thank P.S.Gorain for discussions. 
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